Weak localization, Aharonov-Bohm oscillations and decoherence 

in arrays of quantum dots 
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Combining scattering matrix theory with non-linear cr-model and Keldysh technique we develop a 
unified theoretical approach enabling one to non-perturbatively study the effect of electron-electron 
interactions on weak localization and Aharonov-Bohm oscillations in arbitrary arrays of quantum 
dots. Our model embraces (i) weakly disordered conductors (ii) strongly disordered conductors 
and (iii) metallic quantum dots. In all these cases at T — > the electron decoherence time is 
found to saturate to a finite value determined by the universal formula which agrees quantita- 
tively with numerous experimental results. Our analysis provides overwhelming evidence in favor of 
electron-electron interactions as a universal mechanism for zero temperature electron decoherence 
in disordered conductors. 



I. INTRODUCTION 

Quantum interference of electrons is a fundamen- 
tally important phenomenon which can strongly elec- 
tron transport in disordered conductors^— . Quantum 
coherent effects are mostly pronounced at low tempera- 
tures in which case certain interaction mechanisms are 
"frozen out" and, hence, do not anymore limit the abil- 
ity of electrons to interfere. However, there exists at least 
one mechanism, electron-electron interactions, which re- 
mains important down to lowest temperatures and may 
destroy quantum interference of electrons down to T = 0. 
In a series of papers^ two of the present authors for- 
mulated a general theoretical formalism which allows to 
describe electron interference effects in the presence of 
disorder and electron-electron interactions at any tem- 
perature, including the limit T — > 0. This approach 
extends Chakravarty-Schmid description 2 of weak local- 
ization (WL) and generalizes Feynman- Vernon path in- 
tegral influence functional technique^ to fermionic sys- 
tems with disorder and interactions. With the aid of this 
technique it turned out to be possible to quantitatively 
explain low temperature saturation of WL correction to 
conductance SG WL (T) commonly observed in diffusive 
metallic wires&£. It was demonstrated^ that this satura- 
tion effect is caused by electron-electron interactions. 

It is worth pointing out that low temperature satura- 
tion of WL correction and of the electron decoherence 
time t v (extracted from SG WL (T) or by other means) 
has been repeatedly observed not only in metallic wires 
but also in virtually any type of disordered conductors 
ranging from individual quantum dots^ to very strongly 
disordered 3d structures and granular metals 9 . Hence, it 
is plausible that in all these systems we are dealing with 
the same fundamental effect of electron-electron interac- 
tions. In order to test this conjecture it is necessary to de- 
velop a unified theoretical description which would cover 
essentially all types of disordered conductors. Although 
the approach^ is formally an exact procedure treating 
electron dynamics in the presence of disorder and inter- 
actions, in some cases, e.g., for quantum dots and gran- 



ular metals, it can be rather difficult to directly evaluate 
SG WL (T) within this technique. 

One of the problems in those cases is that the descrip- 
tion in terms of quasiclassical electron trajectories may 
become insufficient, and electron scattering on disorder 
should be treated on more general footing. In addition, 
within the approach- disorder averaging is (can be) post- 
poned until the last stage of the calculation which is con- 
venient in certain physical situations. In other cases - 
like ones studied below - it might be, in contrast, more 
appropriate to perform disorder averaging already in the 
beginning of the whole analysis. Finally, it is desirable to 
deal with the model which would embrace various types 
of conductors with well defined properties both in the 
long and short wavelength limits. 

Below we will elaborate an alternative approach which 
combines the scattering matrix and Keldysh techniques 
with the description of electron-electron interactions in 
terms of quantum Hubbard-Stratonovich fields. Note 
that previously a similar type of approach was employed 
in order to describe Coulomb effects in tunnel junctions, 
see, e.g i 10 ' 11 . Here we will describe a disordered con- 
ductor by means of an array of (metallic) quantum dots 
connected via junctions (scatterers) with an arbitrary 
distribution of transmissions of their conducting chan- 
nels. This model will allow to easily crossover between 
the limits of granular metals and those with point-like 
impurities and to treat spatially restricted and spatially 
extended conductors within the same theoretical frame- 
work. Electron scattering on each such scatterer will 
be treated within the most general scattering matrix 
formalis m 12 ' 13 adopted to include electron-electron in- 
teraction effects^—. Averaging over disorder will be 
performed within the non-linear a— model technique in 
Keldysh formulation. This method has certain advan- 
tages over the imaginary time approach since it allows to 
treat both equilibrium and non-equilibrium problems and 
also enables one to include Coulomb interaction between 
electrons in a straightforward manner—. 

In this paper we will review and extend our analysis of 
weak localization effects and Aharonov-Bohm oscillations 
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in systems composed of metallic quantum dots22r— . In 
Sec. 2 we will construct a theory for essentially non- 
interacting electrons including interaction effects only 
phenomenologically by introducing an effective electron 
dephasing time t v as an independent parameter. In Sec. 
3 we will develop a systematic unified analysis of the ef- 
fect of electron-electron interactions on weak localization 
and Aharonov-Bohm oscillations in both quantum dots 
and extended diffusive conductors. Sec. 4 is devoted to 
a comparison of our results with experimental observa- 
tions. 



II. WEAK LOCALIZATION IN QUANTUM DOT 
ARRAYS 



A. The model and basic formalism 

Let us consider a Id array of connected in series chaotic 
quantum dots (Fig. [I]). Each quantum dot is character- 
ized by its own mean level spacing 5 n . Adjacent quan- 
tum dots are connected via barriers which can scatter 
electrons. Each such scatterer is described by a set of 
transmissions of its conducting channels (here k la- 
bels the channels and n labels the scatterers) . Below we 
will ignore spin-orbit scattering and focus our attention 
on the case of Id arrays. If needed, generalization of our 
analysis to systems of higher dimensions can be employed 
in a straightforward manner—. 

An effective action S[Q] of an array depicted in Fig. 
Q] depends on the fluctuating 4x4 matrix field a 19 i 23 
Qn{t\,t2) defined for each of the dots (n = 1, ...,N — 1). 
Each of these fields is a function of two times t\ and t<x 
and obeys the normalization condition 



1. 



(1) 



The action of an array can be represented as a sum of 
two terms 



iS[Q] = iS d [Q] + iS t [Q]. 



(2) 



The first term, iSd[Q], describes the contribution of bulk 
parts of the dots. This term reads 



N-l 

iS d [Q] = f Tr 



I* 



,H 2 ([A,Q n }Y 



(3) 





FIG. 1: Id array of N— 1 quantum dots coupled by N barriers. 
Each quantum dot is characterized by mean level spacing 5 n ■ 
Each barrier is characterized by a set of transmissions of its 

( n ) 

conducting channels Tj" . 



The second term in Eq. ((2J, iSt[Q], describes electron 
transfer between quantum dots. It has the form 28 



1 N 

i5 «[«] = oEE Trln 



= 1 k 



T, 
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-({Q„-i,Q„}-2) 



.(5) 



Note that here the magnetic field H is included only in 
the term ([3]) describing the quantum dots while it is ig- 
nored in the term ([5]). Usually this approximation re- 
mains applicable at not too low magnetic fields. 

An equilibrium saddle point configuration A(ii — £2) 
of the matrix field Q(t\,t2) depends only on the time 
difference and has the form 



A(t) 



dE 
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(6) 



9 K (E) -1 J 



where g K (E) = 2[1 - 2f F (E)} = 2tanh( J B/2T). This 
choice of the saddle point corresponds to the following 
structure of the 4x4 matrix Green function G: 



G = 



I G A \ 
TG A *T 
—G K G R 
\ TG K *T TG R *T J 



Here we defined the time inversion operator T: 
Tf(t) = f(t f - t), 



(7) 



(8) 



where tf will be specified later. Note that the function 
G in Eq. ([7]), defined for a given disorder configuration, 
should be contrasted from the Green function 



Gn = 



.d_ 

'dt 



2m 



'It, 



Q 



(9) 



Here H is an external magnetic filed, a n = 
b n (e 2 1 'h 2 c 2 )v pd 2 n min{Z e , d„}, b n is a geometry dependent 
numerical prefactor— , d n is the size of n— th dot, l e is the 
elastic mean free path in the dot, and A is 4 x 4 matrix: 



A = 



1000 
0-100 
10 
-1 



(4) 



defined for a given realization of the matrix field Q. In 
Eq. ([9]) we also introduced the electron elastic mean free 
time Tp. 



B. Gaussian approximation 

In order to evaluate the WL correction to conductance 
we will account for quadratic (Gaussian) fluctuations of 
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the matrix field Q n . This approximation is always suffi- 
cient provided the conductance of the whole sample ex- 
ceeds e 2 /h, in certain situations somewhat softer applica- 
bility conditions can be formulated. Expanding in powers 
of such fluctuations we introduce the following parame- 
terization 



Qr. 



= A + i[W n ,A}+W n AW n --{WlA} + O(W 3 ).(10) 

It follows from the normalization condition ([I]) that only 
8 out of 16 matrix elements of W are independent pa- 
rameters. This observation provides certain freedom to 
choose an explicit form of this matrix. A convenient pa- 
rameterization to be used below is 



W n = 




(11) 



With this choice the quadratic part of the action takes 
the form 



iS& = iS2>[a,b]+iS£>[u,v] 



(12) 



(2) 

where iS^ b [a, b] does not depend on H and describes dif- 

fuson modes, while iSuv [u, v] is sensitive to the magnetic 
field and is responsible for the Cooperons. The diffu- 

(2) 

son part of the action iS^J [a, b] was already analyzed 
before^ and will be omitted here. Below we will focus 
our attention on the Cooperon contribution which reads 



n—1 

N-l 



— [U ln ,U 2 r 

at 



n— 1 
N 

^ 2 

n=l 



16a n H u\u 2 



16a n H viv 2 



+ (VlT, 



(Ux n - Ul,„_l)(u 2 
^l,n — 1 )(v2n ~ V2,n-l) 



W2,n-l) 

(13) 



where g n — ^J^fe^fc™'' = ^b,/e 2 R n is the dimensionless 
conductance of n— th barrier. With the aid of the action 
(IT3| we can derive the pair correlators of the fields m >2 
and v\,2- 

(u ln (t 1 ,t 2 )u 2m (t',t")) = (v ln (t',t")v 2m (ti,t 2 )} 

= %±8(tt -h+t'- t")C nm {t" - h), (14) 

where we defined a discrete version of the Cooperon 
C nm (t) obeying the equation 

— 9nC n -l,m — <?n+lCVi+l,m] = S nm S(t). (15) 



This equation should be supplemented by the bound- 
ary condition C nm (t) = which applies whenever one 
of the indices n or in belongs to the lead electrode. Here 
T~Hn = l/l6a n H 2 is the electron dephasing time due to 
the magnetic field. In Eq. (fT5j) we also introduced an ad- 
ditional electron decoherence time in n— th quantum dot 
Ttpn which can remain finite in the presence of interac- 
tions. In this section we account for electron decoherence 
only phenomenologically by keeping the parameter T vn 
in the equation for the Cooperon. Rigorous description 
of quantum decoherence by electron-electron interactions 
will be carried out in Sec. 3. 



C. Weak localization corrections to conductance 

Let us now derive an expression for WL correction to 
the conductance in terms of the fluctuating fields u and v. 
In what follows we will explicitly account for the discrete 
nature of our model and specify the WL correction for a 
single barrier in-between two adjacent quantum dots in 
the array. 

We start, however, from the bulk limit, in which case 
the Kubo formula for the conductivity tensor a a p reads 



<r a p\r,r) = -i l dt' {t-t') 

J — oo 

x (jp(t',r')j a {t,r)-j a (t,r)j (t',r')). (16) 

Following the standard procedure 1,2 , approximating the 
Fermi function as —dfF(E)/dE S(E) (which effec- 
tively implies taking the low temperature limit) and using 
a phenomenological description of interactions as medi- 
ated by external (classical) fluctuating fields 29 , from Eq. 
(|16p one can derive the WL correction in the form: 



x (V^ - V« ) ri=r3=7 .(V^ - v^v i= 

x (G R (t,r 1 ;t"y 2 )G A (t',r' 1 ;t,r 2 )) d 



is, max cross 



(17) 



which implies summation over all maximally crossed di- 
agrams, as indicated in the subscript. At the same time, 
averaging over fluctuations of Q within Gaussian approxi- 
mation is equivalent to summing over all ladder diagrams. 
Since we are not going to go beyond the above approxi- 
mation, we need to convert maximally crossed diagrams 
in Eq. (fT7|) into the ladder ones. Technically this con- 
version can be accomplished by an effective time reversal 
procedure for the advanced Green function which can be 
illustrated as follows. 

Consider, e. g., the second order correction to G A in 
the disorder potential Udi s (x) 



6^G A (t',r' 1 ;t, 



r 2 



i J dr 2 J dri J d x 2 d ,X\ 

x G A (t', r / 1 ;T ll x 1 )U dis (x 1 )G A (T 1 ,x 1 ;T 2 ,x 2 ) 
xU dis (x 2 )G A (T 2 ,x 2 ;t,r 2 ). (18) 
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Making use of the property G A (X 1 ,X 2 ) = G R *(X 2 ,X 1 ), 
we get 

5 {2) G A {t',r' 1 ;t,r 2 )=-i / dr 2 / dn / d 3 x 2 d 3 Xl 
Jt' Jt' J 

x G R *(t, r 2 ;r 2 , a: 2 )lf dis (a; 2 )G K *(T 2 , cc 2 ; n, £1) 
x^diB(xi)G fl *(7i,xi;t',ri). (19) 
Setting i/ = t + 1', we rewrite this expression as follows 

6 (2) G A {t',r[;t,r 2 ) = -i [ ' dr 2 f " dn 
Jt f -t Jtf—t 

x J d 3 x 2 d 3 Xl G R *{t f ~t',r 2 -T 2 ,x 2 ) 



x U dis (x 2 )G R *(T 2 ,x 2 ;n,Xi) 
xU dis {xi)G R *(T 1 ,x 1 ;t f - t,r[). 



(20) 



Close inspection of the right hand side of Eq. (|20p allows 
to establish the following relation 

S^G A (t',r[;t,r 2 ) = T8^G R * (t' ,r 2 ;t,r[)T, (21) 

which turns out to hold in all orders of the perturbation 
theory in Udis- As before, the time inversion operator T 
is defined in Eq. © with t f =t + t'. 

As a result, the expression for Sa^ L takes the form: 

s^(ry) = -^jyjdt" 

x (G R (t, n-X, r 2 )TG R * (f , r 2 ; t, rj )T) dis _ laddcr (22) 

Rewriting Eq. (|22[) in terms of the matrix elements of 
the Green function |[7J), we obtain 



-2 /•* 



dt / di 



x (V° - V« ) ri=ra=P (V^ - V^) K=ri=P / 
x (G33(t,r 1 ;t",r^G44(t',r 2 ;t,r' 1 )) dis laddCT (23) 

Our next step amounts to expressing WL correction 
via the Green function Gq . For that purpose we will 
use the following rule of averaging 

(G33(t,r 1 ;f // ) r- 2 )G 4 4(t',r 2 ;t,r;)> dis 

= (G 33 ., Q (t, n;t", r' 2 )G 4 A;Q(t', r 2 ; t, r[)) Q 

- (G 3 4. Q (t,r 1 ;t,r' 1 )G 43 . Q (t',r 2 ;t",r 2 )) Q . (24) 

One can check that within our Gaussian approximation 
in u and v the first term in the right hand side of Eq. 
(f24f does not give any contribution. Hence, we find 



= 2 r t 



dt' / dt" 



— oo 

x (G 3 4;o(t,r 1 ;t,ri)G43 i Q(* / ,r 2 ;t",ri)) Q . (25) 



Let us now turn to our model of Fig. 1 in which case 
the voltage drops occur only across barriers. In this case 
Eq. (|25|) . which only applies to bulk metals, should be 
generalized accordingly. Consider the conductance of an 
individual barrier determined by the following Kubo for- 
mula 



G 



dt'(t-t'){I(t',x')I(t,x) 



I(t,x)I(t',x')). 



(26) 



Here I(t, x) is the operator of the total current flowing 
in the lead (or dot) and a; is a longitudinal coordinate 
chosen to be in a close vicinity of the barrier. Due to 
the current conservation the conductance G should not 
explicitly depend on x and x'. Comparing Eqs. (|26j) 
and (|16p . and making use of Eq. (|2"5|) and the relation 
I(t,x) = J d 2 zj x (t,x,z), where j x is the current density 
in the direction and z is the vector in the transversal 
direction, we conclude that WL correction to the conduc- 
tance of a barrier between the left and right dots should 
read 



dt' / dt" / d 2 zd 2 z' 



d&LR ~ 4W 

X \^Wl ~~ ^ X2 )xi= X2— X x' ± ~~ ^ x' 2 )x' x — x' 2 —X' 

x {G u -,Q{t, xi,z; t, x' 11 z')G i3 . Q {t', x 2 , z; t", x' 2 ,z')) Q . 

(27) 

In what follows we will assume that both coordinates 
x and x' are on the left side from and very close to the 
corresponding barrier. Let us express the Green function 
in the vicinity of the barrier in the form 

G Q (t, x, z; t', x', z') = {e^-^'G+tAt, t' , x, x') 



+ e- i ^+ i ^'0--(t,t',x,x') 
+ e i ^+ i ^*'$+-(t,t',x,x') 



-ip n x — ip m x 



g-+(t,t',x,x')}<Z> n (z)$* m (z'), (28) 



where ^ n {z) are the transverse quantization modes 
which define conducting channels, p n is projection of 
the Fermi momentum perpendicular to the surface of the 
barrier, and the semiclassical Green function Q^ n slowly 
varies in space. Eq. (|27|) then becomes 



6GY R L 



dt' / dt" 



E 



E 

mnkl a/37<5— ±1 



(ap n - jp k )(/3p m - Spi) 



x e 



(G^l M (t,t,x,x')gif. A3 (t',t'',x ) x')) i 

iap n xi—i0p rn x 1 -\-i-ypkX2—'iSpi x' 2 



(29) 



Next we require SG^^ to be independent on x and 
x' , i.e. in Eq. (|29[) we omit those terms, which contain 
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quickly oscillating functions of these coordinates. This 
requirement implies that ap n + "/p^ = and f3p m + 5pi = 
0. These constraints in turn yield 7 = —a, i5 = — /3, 
k = n and ^ = m. Thus, we get 



-,1 



E 



mn a/9=±l 



di' / dt" CtfipnPn 



X (^;34(*.*^^0^:4t(*'>*''^^')) Q - (so) 

Let us choose the basis in which transmission and 
reflection matrices t and r are diagonal. In this ba- 
sis the semiclassical Green function is diagonal as well, 
Qmn oc GnnSnm, and Eq. f[30]) takes the form 



P 2 . r> 2 /"* 

-E^/ 

„ ™ ./-oc 



L,nn:34 



L,nn;43 



(*',*") 

^L,nn;34(^ *)^L,nn;43(* ' ^ ) 
^7 m ;34(M)eZ+„;4 3 (^0 
L.nn;34 ,nn;43 



-g- + 



Q' 



(31) 



What remains is to express the WL correction in terms 
of the field Q only. This goal is achieved with the aid of 
the following general relation^ 3 - 



dt' / dt" 



x (T n [v 1L (t,t)v 2R (t',t") + v 1R (t,t)v 2L (t',t")] 

+ T%[v 1L (t, t) - v 1R (t, t)} [v 2L (f, t") - v 2R (t', t")}). (32) 

Note that the contribution linear in T n , which contains 
the product of the fluctuating fields on two different sides 
of the barrier, vanishes identically provided fluctuations 
on one side tend to zero, e.g. if the barrier is directly 
attached to a large metallic lead. In contrast, the contri- 
bution oc T 2 in Eq. (|32|) survives even in this case. 
Finally, applying the contraction rule (fT4")) we get 

p 2 poo 

SG^r = J dt{p[5 R C LR (t) + 5 L C RL (t)] 

+ (1 - (3)[S R C RR (t) + 8 L C LL {t)\ }. (33) 

Here 8l jR is the mean level spacing in the left/right quan- 
tum dot, 
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(34) 



is the dimensionless conductance of the barrier and 

P = J2Tk(l-T k )/J2n (35) 

k k 

is the corresponding Fano factor. 

Likewise, the WL correction to the n— th barrier con- 
ductance in Id array of N — 1 quantum dots with mean 




FIG. 2: Single quantum dot connected to the leads via two 
barriers. 



level spacings 8 n connected by N barriers with dimen- 
sionless conductances g n and Fano factors j3 n reads 



8G 



WL 



' Air 2 



dt{/3 n [8 n C n -i tn (t) 



+ «n-lCn,T.-l(t)] + (1 - Pn) [S n C nn {t) 

+ <5»-iC»-i,„-i (*)]}. (36) 

So far we discussed the local properties, namely WL 
corrections to the conductivity tensor, 8cr^g(r,r'), and 
to the conductance of a single barrier, 8G^^ . Our main 
goal is, however, to evaluate the WL correction to the 
conductance of the whole system. For bulk metals one 
finds that at large scales the WL correction (fl7|) is lo- 
cal, 6a^p(r,r') oc 8(r — r'). In general though, there 
can exist other, non-local, contributions to the conduc- 
tivity tensor—. Without going into details here, we only 
point out that, even if these non-local terms are present, 
one can still apply the standard Ohm's law arguments 
in order to obtain the conductance of the whole sam- 
ple. Specifically, in the case of Id arrays one finds 23 (sec 
also2i) 



8G 



WL 



1 



1 



sg^y 



En=^Gr/9i 



higher order terms. (37) 



Eqs. J33J), (|36J) and §7$ will be used to evaluate WL 
corrections for different configurations of quantum dots 
considered below. 



D. Examples 

1. Single quantum dot 

We start from the simplest case of a single quantum 
dot depicted in Fig. 2. In this case the solution of Eq. 
(fIBl reads 



Cn(t) = exp 



t 

TD 



t 
TH 



t 



(38) 



where t d — iir/(gi + g 2 )8^ is the dwell time, and 8 c i is 
the mean level spacing in the quantum dot. All other 
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FIG. 3: Most general system with two quantum dots 



components of the Cooperon are equal to zero. From 
Eq. J33J we get 



8G{ 



WL 



e 2 g 1 (l-j3 1 )S d 



1 



SG 



WL 



4tt 2 1/td + 1/t h + 1/t v ' 
e 2 g 2 (l-f3 2 )5 d 1 



.(39) 



4tt 2 1/td + 1/t h + 1/t v 
According to Eq. (|3T|) the total WL correction becomes 

wl e 2 S 5iff 2 2 (l-/3i)+.9?.92(l-/3 2 ) 



SG 



4tt 2 ( gi + g 2 ) 2 (1/t d + 1/t v + 1/th) ' 



(40) 



Since \/t h oc H 2 , the magnetoconductance has the 
Lorentzian shaped. In the limit H = and in the ab- 
sence of interactions (t v —> 00) Eq. (|4"0")) reduces to^ 2 . 



5G 



WL 



e 2 g l g 2 {l-p l )+g 2 g 2 {l-p 2 ) 
* (91+92) 3 



(41) 



As one can see for the case of low transmissions (for ex- 
ample in case of tunneling barriers) the WL corrections 
equals to zero. 



2. Two quantum dots 

Next we consider the most general setup composed of 
two quantum dots with the corresponding conductances 
and Fano factors defined as in Fig. [3l The Cooperon 
is represented as a 2 x 2 matrix which zero frequency 
component satisfies the following equation 



5ii + .912 + g y + 7i " 9y \ ( G11 C12 

~9y 921 + 522 + g y + 72 / V ^ 22 

( ^/Sl 






A-k/8-2 



where 



7l,2 



47T 

£1,2 



Tfll,2 



Twl, 2 



(42) 



(43) 



Defining A = (3ii+3i2+.9y + 7i)(32i+.922+.9y + 72)-.9y, 
we get 



(521 +322 + 9 V +l2)/Si g y /5 2 
9y/5i {gn + 312 + g y + ii)/8 2 





Ci 2 \ 




{ C21 


G 22 1 






FIG. 4: The magnetoconductance of two dots of Fig. [3] for 
di,d 2 » l e , di/d,2 = 5, Qij = go, Pa = 0, f} y = 0, r v \ = 
t V 2 = 00. Here Hi — 1/4-^/aiTm is the field at which weak 
localization is effectively suppressed in the first dot. For g y = 
the magnetoconductance is given by superposition of two 
Lorentzians with different widths (decoupled dots), while for 
large g y only one Lorentzian survives corresponding to the 
contribution of a one "composite dot" . 




FIG. 5: Two quantum dots in series. 



With the aid of Eq. (|33p we can derive WL corrections 
for all five barriers in our setup which we do not specify 
here for the sake of brevity (see 2 ^ for further details). 

WL correction to the conductance of the whole struc- 
ture 8G WL is obtained from the general expression for 
the conductance determined by Ohm's law: 

G = [GnGi 2 (G 2 i + G 22 ) + G 2 iG 22 (Gn + Gi 2 ) 
+ G y (Gi2 + G 22 )(Gn + G21)] 
/ [(Gn + Gi 2 )(G 2 i + G 22 ) 
+ G y (G n + G 12 + G 21 + G 22 )). (44) 

Substituting Gy — > Gy + SG^ L into this formula and 
expanding the result to the first order in SG^ L , we get 



sg wl = J2 J^LsG™ L + ^5G™ L 



i,J=l,2 



dG t 



dG y v 



(45) 



This general result for the WL correction to the conduc- 
tance is illustrated in Fig. [4] for a particular choice of the 
system parameters. 

Of particular importance for us here is the system of 
two quantum dots connected in series, as shown in Fig. [SJ 
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i.e. in the general structure of Fig. [3]we set G12 = G21 = 

0, Gn — Gl, Gy = G2, G22 = G3, PH = Pl, Py = /?2 

and P22 = p3- We also assume H = and = 00. 
WL corrections to the barrier conductances then take 
the form 



WL 



5G{ 



8Gf L 



5Gf L 



e 2 gi(g2 + g 3 )(l-/3i) 

7T 5152 + 5253 + 5153 ' 

e 2 ff2(gi+g 3 )(l-/? 2 ) + 2gl 

7T glg2 + 5253 + glg3 

e 2 g3(gi +52XI- As) 



5152 + 3233 + 51.93 

while Eq. (|4"4")l reduces to 



G 



G1G2G3 



G1G2 + G1G3 + G2G3 



(46) 



(47) 



WL correction for the whole system then reads 



<5G 



WL 



e 2 5i52 2 5l(52 +53)(l-fii) 

7T (5152 + 5253 + 5153) 3 

e 2 5i5 2 5l(5i+53)(l-/3 2 ) 
7T (5152 + 5253 + 5i 53) 3 
e 2 5 2 5 2 2 53(5i+52)(l-/3 3 ) 



7T (5l52 

2e 2 



-5253 + 5i 53 J 

2 2 2 

5i525s 



7f (5152 + 5253 + 5i53) 3 



(48) 



In the limit of open quantum dots, i.e. Pi, 2, 3 = 0, we 
reproduce the result^. It is easy to see that provided 
the conductance of one of the barriers strongly exceeds 
two others, Eq. (1451 reduces to Eq. fill . If all three 
barriers are tunnel junctions, /3i,2,3 —> 1, the first three 
contributions in Eq. f|48jl vanish, and only the last con- 
tribution - independent of the Fano factors - survives in 
this limit. If, on top of that, one of the tunnel junctions, 
e.g. the central one, is less transparent than two oth- 
ers, <?2 <C 51,53, the result acquires a particularly simple 
(non-Lorentzian) form 



SG WL = _ 



2e 2 



53 



^ (51 +7l) (53 +72)' 



(49) 



with 71,2 defined in Eq. fig]). Note that SG WL oc g%, i.e. 
this result is dominated by the second order tunneling 
processes across the second barrier. 



and T vn = Tip. For this system the Cooperon can also be 
found exactly. The result reads 

N-i 



C nm (u)) 



N ^ 

9=1 



TTOn • Tram 

sin -ft- sin -f r 



TH Tu, 



(50) 



Here td = 2-KjgSd and th = 1/lQaH 2 . The WL correc- 
tion then takes the form 



5Q WL = 



2 <r N ~ x 

e gSd 



(3 cos 22 



N 



2n 2 N 2 ^ J_ 

9=1 TH 



(51) 



The sum over q can be handled exactly and yields 



SG 



WL 



ttN 2 



N- 



.1 



,27V 



1 



1 



,2N 



1-U 2 



13(1 + u 2 ) + 2(1 - p)u 
1-u 2 



(N-l)/3 



where 



u = l 




(52) 



(53) 



In the tunneling limit (3 = 1 and for t v — > 00 our result 
defined in Eqs. (|52 |) -([53 | becomes similar - though not 
exactly identical - to the corresponding result 33 . 

If Tip is long enough, namely 1/t v < i/rh, where -&rh = 
7r 2 /2N 2 td is the Thouless energy of the whole array in 
Eqs. (|5"T j) - (|5li|) it is sufficient to set t v — 00. In this 
case the magnetic field H significantly suppresses WL 
correction provided 1/rjj > £rh or, equivalently, if 



H > H 



N • 



H 



N 



1 

8N 



irg5 d 



(54) 



In the opposite limit 1 /t v > £xh we find 



SG WL = 



ttN 



/?(l + ^ + ^ 



l + ln + ln 

TH Tu, 



.(55) 



In particular, in the diffusive limit th,t v 3> td we get 



SG 



WL 



Dt h t v 



TlNd 



(56) 



3. ID array of identical quantum dots 

Let us now turn to Id arrays of quantum dots depicted 
in Fig. [TJ For simplicity, we will assume that our array 
consists of N — 1 identical quantum dots with the same 
level spacing 5 n = Sd and of N identical barriers with the 
same dimensionless conductance g n = g and the same 
Fano factor (3 n = j3. We will also assume that the quan- 
tum dots have the same shape and size so that TH n = t h 



where we introduced the diffusion coefficient 

D = d 2 /2T D . (57) 



Eq. (|56|) coincides with the standard result for quasi- 
ld diffusive metallic wire. Note, however, that the val- 
ues of th within our model may differ from those for a 
metallic wire. The ratio of the former to the latter is 
T ij d / r H Ct ~ T fi/ T D, where ta ~ d/vp is the flight time 
through the quantum dot. Since typically ta < td we 
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FIG. 6: Magnetoconductance of a Id array of N — 1 identical 
open (/3 = 0) quantum dots in the absence of interactions 
[t v — > oo). The field Hn is defined in Eq. (53) . 



conclude that for the same value of D the magnetic field 
dephases electrons stronger in the case of an array of 
quantum dots. 

For a single quantum dot [N = 2) Eq. (j5"2")) reduces to 



SG WL = _ 



e 2 (l-/3) 



47T 




- LB. 






TH 





(58) 



in agreement with Eq. (|40|) . 

For two identical quantum dots in series we obtain 



SG WL = - — 
9?r 



2-/3 



1 



2t d 
th 



2t d 



1 



2r o | 2-rp 
3rff 3t> 



,(59) 



i.e. the magnetoconductance is just the sum of two 
Lorentzians in this case. 

Finally, in the absence of any interactions (t v — oo) 
and at H = we obtain 



SG 



WL 



e 



1/9 1 fa 1 
3~N + N2V~3 



(60) 



In the limit N — > oo this result reduces to the standard 
one for a long quasi- Id diffusive wire^s while for any finite 
N we reproduce the results for tunnel barriers^ (/3 —> 1) 
and open quantum dots^ (/3 — > 0). 

The magnetoconductance of a Id array of N — 1 iden- 
tical open quantum dots in the absence of interactions is 
also illustrated in Fig. 6. 



III. QUANTUM DECOHERENCE BY 
ELECTRON-ELECTRON INTERACTIONS 

A. Qualitative arguments 

Let us now include electron-electron interactions and 
analyze their impact on loss of phase coherence of elec- 




FIG. 7: Single quantum dot and a pair of time-reversed elec- 
tron paths. Fluctuating voltages Vl and Vr are assumed to 
drop only across the barriers and not inside the dot. 



trons' wave functions. Before turning to a detailed calcu- 
lation it is instructive to discuss a simple qualitative pic- 
ture demonstrating under which conditions decoherence 
by electron-electron interactions is expected to occur. 

Consider first the simplest system of two scatterers 
separated by a cavity (quantum dot, Fig. 7) The WL 
correction to conductance of a disordered system Gw l is 
known to arise from interference of pairs of time-reversed 
electron paths 2 . In the absence of interactions for a sin- 
gle quantum dot of Fig. 7 this correction was evaluated 
in the previous sections (see Eq. (HU). The effect of 
electron-electron interactions can be described in terms 
of fluctuating voltages. Let us assume that the voltage 
can drop only across the barriers and consider two time- 
reversed electron paths which cross the left barrier (with 
fluctuating voltage Vt(t)) twice at times U and tf, as it 
is shown in Fig. 7. It is easy to see that the voltage- 
dependent random phase factor exp(i J. f VL(t)dt) ac- 
quired by the electron wave function >F along any path 
turns out to be exactly the same as that for its time- 
reversed counterpart. Hence, in the product <F\]/* these 
random phases cancel each other and quantum coher- 
ence of electrons remains fully preserved. This implies 
that for the system of Fig. 9 fluctuating voltages (which 
can mediate electron-electron interactions) do not cause 
any dephasing. 

This qualitative conclusion can be verified by means 
of more rigorous considerations. For instance, it was 
demonstrated^ that the scattering matrix of the system 
remains unitary in the presence of electron-electron in- 
teractions, which implies that the only effect of such in- 
teractions is transmission renormalization but not elec- 
tron decoherence. A similar conclusion was reached^ 
by directly evaluating the WL correction to the system 
conductance. Thus, for the system of two scatterers of 
Fig. 7 electron-electron interactions can only yield en- 
ergy dependent (logarithmic at sufficiently low energies) 
renormalization of the dot channel transmission a 18 ' 20 but 
not electron dephasing. 

Let us now add one more scatterer and consider the 
system of two quantum dots depicted in Fig. 8. We again 
assume that fluctuating voltages are concentrated at the 
barriers and not inside the cavities. The phase factor 
accumulated along the path (see Fig. 8) which crosses the 
central barrier twice (at times ti and t > U) and returns 
to the initial point (at a time tf) is e l ^ p ^ ti ^ v ^\ where 
(p/e = V(t) is the fluctuating voltage across the central 



FIG. 8: Two quantum dots and a typical electron path. Fluc- 
tuating voltages Vl , V and Vr are again assumed to drop only 
across the barriers. 



barrier. Similarly, the phase factor picked up along the 
time-reversed path reads e*[ ¥> (* J ' +t< ~ t ) _ ' p Wl. Hence, the 
overall phase factor acquired by the product ^<F* for a 
pair of time- reversed paths is exp(z<I>tot), where 

*tot (*<,*/,*) = <p(ti) ~ <p(t) - <f + {t f +t 4 -t) + <p(tf). 

Averaging over phase fluctuations, which for simplicity 
are assumed Gaussian, we obtain 

^ e i<E>*o*(fi>*/>*)\ = e - 1 <*?<,*(*< .*/.*)} 

_ e -2F(t-U)-2F(t f -t)+F(t f -U)+F(t f +U-2t)^ /g-g 

where we defined the phase correlation function 

F(i) = ((^)-^(0)) 2 )/2. (62) 

Should this function grow with time the electron phase 
coherence decays and, hence, Gwl has to be suppressed 
below its non-interacting value due to interaction- 
induced electron decoherence. 

The above arguments are, of course, not specific to sys- 
tems with three barriers only. They can also be applied to 
any system with larger number of scatterers, i.e. virtually 
to any disordered conductor where - exactly for the same 
reasons - one also expects non-vanishing interaction- 
induced electron decoherence at any temperature includ- 
ing T = 0. Below we will develop a quantitative theory 
which will confirm and extend our qualitative physical 
picture. We are going to give a complete quantum me- 
chanical analysis of the problem which fully accounts for 
Fermi statistics of electrons and treats electron-electron 
interactions in terms quantum fields produced internally 
by fluctuating electrons. 



B. Nanorings with two quantum dots 

1. The model and basic formalism 

Consider the system depicted in Fig. 9. The structure 
consists of two chaotic quantum dots (L and R) charac- 
terized by mean level spacing Sl and S R which are the 
lowest energy parameters in our problem. These (metal- 
lic) dots are interconnected via two tunnel junctions Ji 
and J2 with conductances Gti and Gt2 forming a ring- 
shaped configuration as shown in Fig. 9. The left and 



FIG. 9: Two quantum dots with magnetic flux. 

right dots are also connected to the leads (LL and RL) re- 
spectively via the barriers Jl and 3 R with conductances 
Gl and G R . We also define the corresponding dimension- 
less conductances of all four barriers as <?ti,2 = GtipRq 
and gL.R — Gtx.iRq, where R q = 2ir/e 2 is the quantum 
resistance unit. 

The whole structure is pierced by the magnetic flux $ 
through the hole between two central barriers in such way 
that electrons passing from left to right through different 
junctions acquire different geometric phases. Applying a 
voltage across the system one induces the current which 
shows AB oscillations with changing the external flux 
Note that in the absence of the magnetic flux the system 
just reduces to that of two connected in series quantum 
dots (cf. Fig. 5) which is also subject to weak localiza- 
tion effects. Thus, the model considered here allows to 
analyze WL and AB effects within the same formalism 
to be developed below. The system depicted in Fig. 9 is 
described by the effective Hamiltonian: 

H = ^2 l} 2 3 + Hll + Hrl 

i,j=L,R 

+ J2 H 3 +f L +f R + f, (63) 

j=L,R 

where Cij is the capacitance matrix, Vl(r) 1S the electric 
potential operator on the left (right) quantum dot, 

Hll= ]T [ Ml LL (r)(H LL -eV LL )i> a , LL (r), 

H RL =Y. I d 3 r%l tRL (r){H RL -eV RL )V a< R L (r) 

are the Hamiltonians of the left and right leads, Vll.rl 
are the electric potentials of the leads fixed by the exter- 
nal voltage source, 




defines the Hamiltonians of the left (j — L) and right 
(j = R) quantum dots and 
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is the one-particle Hamiltonian of electron in j-th quan- 
tum dot with disorder potential Uj (r) . Electron transfer 
between the left and the right quantum dots will be de- 
scribed by the Hamiltonian 

£ / <*V[t(r)$^(r)$ a , fl (r) + c.c.]. 



a=t4' 



The Hamiltonian describing electron transfer be- 

tween the left dot and the left lead (the right dot and the 
right lead) is defined analogously. 

The real time evolution of the density matrix of our 
system is described by means of the standard equation 



Pit) 



e~ iHt p e* Ht , 



(64) 



where H is given by Eq. (|63[) . Let us express the oper- 
ators e~ tHt and e lHt via path integrals over the fluctu- 



F B 

ating electric potentials Vj ' defined respectively on the 
forward and backward parts of the Keldysh contour: 

- f DVf T exp i-i I" dt'H [vf(t')] ) , 



< = / DVf T exp \ i J dt'H [vf(t')] } . (65) 



Here T exp (T exp) stands for the time ordered (anti- 
ordered) exponent. 

Let us define the effective action of our system 



iS[V F ,V B ] = In tr 



T exp<^ -i / dt'H [Vfit')] 



/A i T c.xp<J i J dt'H [V 3 B (t')\ 



66) 



Integrating out the fermionic variables we rewrite the 
action in the form 



iS = iS c + iS ext + 2Trln [GT 1 ] 



(67) 



Here Sc is the standard term describing charging effects, 
Sext accounts for an external circuit and 



(Gil T L 



G" 1 



V o 





rjf &1 1 f o 

ft Gn 1 f R 

R 







(68) 



T " G RL f 



is the inverse Green-Keldysh function of electrons prop- 
agating in the fluctuating fields. Here each quantum dot 
as well as two leads is represented by the 2x2 matrix in 
the Keldysh space: 



G7 l = 



id. 



Hi 




eV/ 





-id t + Hi - eV/ 



(69) 



Ct3 





FIG. 10: Diagrammatic representation of different contribu- 
tions originating from expansion of the effective action in pow- 
ers of the central barrier transmissions: second order (AES) 
terms (a) and different fourth order terms (b,c). 



2. Effective action 

Let us expand the exact action iS (|67[) in powers of T. 
Keeping the terms up to the fourth order in the tunneling 
amplitude, we obtain 



iS w iSc + iS e 



iS L + iS R - 2tr 



GrfGRl* 



-tr 



G L fG R f^G L fG R f^ 



(70) 



Here iS l jR are the contributions of isolated dots, the 
terms oc t 2 yield the Ambegaokar-Eckern-Schon (AES) 
action-^ iS AES described by the diagram in Fig. 10a, 
and the fourth order terms oc t (diagrams in Fig. 10b, c) 
account for the weak localization correction to the system 
conductance^^. 

It is easy to demonstrate^ that after disorder averag- 
ing iS AES becomes independent of $ and, hence, it does 
not account for the AB effect investigated here. After 
averaging the last term in Eq. ([70| over realizations of 
transmission amplitudes and over disorder only the con- 
tribution generated by the diagram (c) keeps depending 
on the magnetic flux and yields^ 



mi9t2 



A^N L N R ^ 

m,n— 1.2 



J dridr 2 J dti...dt4,C L (Ti)C R (T2) 



i{ 9 +{t 2 )- 9 +{U)+ip+{U)->p+{H)) ei n t 



xe 



h(h -t 2 - r^e 1 ^—- 



h(t 2 



T 2 )e 



-f(h-t 2 



n)e 



-f{h - ta - r 2 )e l ^ i /i(i 3 -U + n) 



. ip~ (t d ) 



■ y (t 4 ) 

-e h(U-t l +T 2 ) 



-{L o R, ip d 



(71) 
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where CL,R,(t) the Cooperons in the left and right dots, 
f(t)= j f F (E)dE/2TT is the Fourier transform of the 
Fermi function Jf{E) and h(t) = 8(t) — f(t). Here we 
also introduced the geometric phases 



R 

^ J dx^A^x), 

L 



(72) 



where the integration contour starts in the left dot, 
crosses the first {pp) or the second {<Pp) junction and 
ends in the right dot. The difference between these two 
geometric phases is ipg — ipP — 27r$/<I>o. In addition, 
we defined the "classical" and the "quantum" compo- 
nents of the fluctuating phase <p + (t) = {(pp(t) + <pB(t))/2 
and <f~{t) = ^PF{t) — <fB(t), where the phases (Pf,b(£) = 
e J dr(Vp' B (r) — Vp B (t)) are defined on the forward 
and backward parts of the Keldysh contour. 

The above expression for the action S^ L (I7T1) fully 
accounts for coherent oscillations of the system conduc- 
tance in the lowest non- vanishing order in tunneling. The 
WL contribution to action of two quantum dots is recov- 
ered in exactly the same way2i The result is the similar 
except geometric phases should be omitted and the com- 
bination gt\Qt2 should be substituted by g 2 x or g\ 2 . 



3. Aharonov-Bohm conductance and WL correction 

Let us now evaluate the current I through our system. 
This current can be split into two parts, I = la + 51, 
where Iq is the flux-independent contribution and 51 is 
the quantum correction to the current sensitive to the 
magnetic flux $. This correction is determined by the 
action iS^ L , i.e. 



51 = 



XQWL! + 

2>V * [ : 



5<p-(t) 



(73) 



In order to evaluate the path integral over the phases 
( ( 9 ± in (1731) we restrict our consideration to the most inter- 
esting for us metallic limit assuming that dimensionless 
conductances gi,,R are much larger than unity, while the 
conductances gn and g t 2 are small as compared to those 
of the outer barriers, i.e. 



9l,9r > 1,5*1,3*2- 



(74) 



In the limit (|74j) phase fluctuations can be considered 
small down to exponentially low energie s 14 ' 37 in which 
case it suffices to expand both contributions up to the 
second order <^ ± . Moreover, this Gaussian approximation 
becomes exac t 15 ' 18 ' 20 ' 21 in the limit of fully open left and 
right barriers with gh,R 3> 1- Thus, in the metallic limit 
(|74[) the integral (1731 remains Gaussian at all relevant 
energies and can easily be performed. 

This task can be accomplished with the aid of the fol- 
lowing correlation functions 



(<p+{t))=eVt, for(t)) =0, 



(75) 



((<p + {t)-<p + (fl))<p + (Q))=-F(t), 



(76) 



<^+(%-(0) + <p-(t)<p+(0)) = 2iK(\t\), (77) 



(¥> + («)¥>"(0) - ^(% + (0)) = 2iK(t), (78) 



(¥>-(*)¥>-(())) =0, 



(79) 



where the last relation follows directly from the causality 
principle- Here and below we define V = Vrl — Vll to 
be the transport voltage across our system. 

Note that the above correlation functions are well 
familiar from the so-called P(_E)-theor} r 10 ' 38 describing 
electron tunneling in the presence of an external environ- 
ment which can also mimic electron-electron interactions 
in metallic conductors. They are expressed in terms of an 
effective impedance Z(lu) "seen" by the central barriers 
Ji and J2 



F(t) = e 2 J ^coth^[ZM] 



1 — cos(ujt) 



ui 



' Z7T UJ 



(80) 



(81) 



Further evaluation of these correlation functions for our 
system is straightforward and yields 



F(t) ~ - In 
9 



sinh(TrTt) 



2tt 



+ 7 



K(t) ~ — sign(t), 



(82) 



(83) 



where we defined g = 4ir/e 2 Z(0) and 7 ~ 0.577 is the 
Euler constant. Neglecting the contribution of external 
leads and making use of the inequality (|74l) we obtain g ~ 
^■9l9r/ (gz +9r)- We observe that while F(t) grows with 
time at any temperature including T = 0, the function 
K{t) always remains small and it can be safely ignored 
in the leading order in 1/g <C 1. After that the Fermi 
function /f(E) drops out from the final expression for 
the quantum correction to the current^—. Hence, the 
amplitude of AB oscillations is affected by the electron- 
electron interaction only via the correlation functions for 
the "classical" component of the Hubbard-Stratonovich 
phase ip + . 

The expression for the current takes the form 



«($) 



-Iab cos(4tt$/$o) - Iwli - IwL2, (84) 



where the first - flux dependent - term in the right-hand 
side explicitly accounts for AB oscillations, while the 
terms Iwli. 2 represent the remaining part of the quan- 
tum correction to the current^! which does not depend 
on 3>. 
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Let us restrict our attention to the case of two identical 
quantum dots with volume V, dwell time Tp and dimen- 
sionless conductances g^ = gR = g = Att/Std, where 
S = 1/Vv is the dot mean level spacing and v is the elec- 
tron density of states. In this case the Cooperons take 
the form C L (t;x,y) = C R (t;x,y) = (0(i)/V)e"*/ TD . We 
obtain^ 



in the case of diffusive metals by means of a different 
approach 4 . Thus, in order to evaluate the decoherence 
time for interacting electrons in arrays of quantum dots 
it is sufficient to account for the fluctuating fields V + 
totally ignoring the fields V~ . The corresponding calcu- 
lation is presented below. 



l AB 



e 2 g t igt2S 2 V 

47T 3 



d,Tid,T2e T o 



(85) 



I 



wh\a 



8tt 3 



dnd^e- 1 ^ 2 -^ 1 ^ . (86) 



where T = 2F(n) + 2F{t 2 ) - F(n - r 2 ) - F(n + t 2 ). 

In the absence of electron-electron interactions this for- 
mula yields 7^ = 4e 2 gtigt2V/ (jtg 2 ). In order to account 
for the effect of interactions we substitute Eq. (|82|) into 
Eq. (|85l) . Performing time integrations at high enough 
temperatures we obtain 



Iab 

AB 



1 

2t d 



(2ttTt rc ) 



r D l <T<r 



1 

AO 



' RC 



(87) 



while in the low temperature limit we find 



Iab 

AO) 

AB 



TD 



s/g 



T < 



D 



Essentially the same results follow for Iwli.2- These re- 
sults demonstrate that interaction-induced suppression of 
both AB oscillations and WL corrections in metallic dots 
with trc t d persists down to T = 0. The fundamen- 
tal reason behind this suppression is that the interaction 
of an electron with an effective environment (produced by 
other electrons) effectively breaks down the time-reversal 
symmetry and, hence, causes both dissipation and de- 
phasing for interacting electrons down to T = 0^. In 
this respect it is also important to point out a deep re- 
lation between interaction-induced electron decoherence 
and the P(.E)-theor y lcl i 38 which was already emphasized 
elsewhere 24-26 . 



Arrays of quantum dots and diffusive 
conductors 



One of the main conclusions reached above is that the 
electron decoherence time is fully determined by fluctua- 
tions of the phase fields tp + (and the correlation function 
F(t)), whereas the phases <p~ (and the response function 
K(t)) are irrelevant for t v causing only a weak Coulomb 
correction to Gwl- This conclusion is general being in- 
dependent of a number of scatterers in our system. Note 
that exactly the same conclusion was already reached 



1. Id structures 

Let us consider a ID array of N — 1 quantum dots by 
N identical barriers as shown in Fig. 1. For simplicity, 
we will stick to the case of identical barriers (with di- 
mensionless conductance j > 1 and Fano factor j3) and 
identical quantum dots (with mean level spacing 5 and 
dwell time td = 2ir/Sg). The WL correction to the sys- 
tem conductance has the form (see Eq. (|5o) ): 



G 



WL 



e 2 gS 
' Att 2 N 2 



"I JO 



x {P[C n - ltn (t) + C ntn -i(t)\ 

+ (1 - P)[C m (t) + C n - 1>n -i(t)] }. (89) 

The Cooperon C nm (t) is determined from a discrete ver- 
sion of the diffusion equation. For non-interacting elec- 
trons and in the absence of the magnetic field this equa- 
tion reads 



dt 



2t d 



SnmS(t). (90) 



The boundary conditions for this equation are C nm = 
as long as the index n or m belongs to one of the bulk 
electrode. The solution of Eq. (jHOJ) with these boundary 
conditions can easily be obtained. We have 



(7(0) m = A V ( — e~ iut 
nm[ ) N ^ J 27T 

q=l J 



•nan • nam 

sin -ft- sin - J j r - 

, i-cos^a 



(91) 



This solution can be represented in the form Cnm{t) = 
Ci-Jl(t)-C^Jt), where 



•7=1 



cos 



Trq(n-rn) 
N 



1 — cos - 



(92) 



In the limit of large N the term C^ 1 " 1 ^ (t) can be safely 
ignored and we obtain C nm {t) w C^^it). Let us express 
the contribution C^^(t) as a sum over the integer val- 
ued paths v{t), which start in the m— th dot and end in 
the n— th one (i.e. v(0) = m, v(t) = n) jumping from one 
dot to another at times tj . This expression can be recov- 
ered if one expands Eq. (|92|) in powers of t^ 1 cos[jrq/N] 
with subsequent summation over q in every order of this 
expansion. Including additional phase factors acquired 
by electrons in the presence of the fluctuating fields , 
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we obtain 



c nm {t) — 

k=\n—m\ i^(t) 
I rt rth 



i/(t)—n 
i/(0)=m 



, , dt k I dtk-i ■■ dt 2 dti 
(2n>) fc J J a J Jo 



t2—t2 tL 2-tl tl 

x e t d e T o . . . e t -d e T o e T o 



x exp 



Averaging over Gaussian fluctuations of voltages V + 
and utilizing the symmetry of the voltage correlator 
(V+(n)V+(r 2 )) = (V+(n)V+(T 2 )), we get 



«/(*)= 

"(0)= 



c nm (t) — ^2 

fc— |n— m| f("r) 

, — t/TD ft ptk 

dt k I dtk-i 



t3 pt2 

dt 2 / dti 
o Jo 



x exp | - e 2 ^ dn J* dr 2 [(V+^tJV*^ (r 2 )) 

-(v: {Tl) (n)K%_ T2) (r 2 ))]}. (94) 



The correlator of voltages can be derived with the aid 
of the er-model approach developed in Sec. 2 of this pa- 
per. Integrating over Gaussian fluctuations of the Q- 
fields one arrives at the quadratic action for the fluctu- 
ating fields V + which has the form 



._ 2i \ f du 

9=1 J 



AC ( 1 — cos — — - 



C a 



gr D e' 



7T —iLUTD + 1 — COS 5? 



N , 

V+(lj)V-(-u>) 



N 



-Y 



duj gr^e 2 (l - cos $) w coth ^ 



N U J 2n 71 ^ T |, + (l-cosf) 



x V7(w)V-(-a;). 
Here we defined 



(95) 



(96) 



The action (j9"5)) determines the expressions for both cor- 
relators (V+V+) (F-function) and (V + V~) (^-function) 
responsible respectively for decohcrcncc and Coulomb 
blockade correction to WL. Since our aim is to describe 
electron decoherence, only the first out of these two cor- 



relation functions is of importance for us here. It reads 

N-l 



q=l 



du) 
2^ 



^(l-cos$) sin^isin^fi 



AC 1 



•OOS$) +C g 



'JT D < - 



— cos -S 



Tp W COth 



«a^ + (l-co8^ 



(97) 



In the continuous limit N 3> 1 and for sufficiently low 
frequencies a; <C 1/td both correlators and 
(V + V~) defined by Eq. (|95|) reduce to those of a dif- 
fusive metal^. 

To proceed let us consider diffusive paths v(t), in 
which case one has 



1 



N-l 



N-l ^ 

n,m— 1 

(In - t 2 |), 



(^n + (n)K+(r 2 )) 



(98) 

where D nm (r) is the diffuson. For ff — »• it exactly 
coincides with the Cooperon for non-interacting electrons 

(0) 



9=1 J 



7rqn • Tram 

■fa* sm sm 



Substituting Eq. (|98) into (|9l)l. we obtain 



C nTO (*)«CW(t). 



-jf(*) 



(99) 



(100) 



where 

J-(t) = 



N-l 
\ " 



N-l ^ 



n.m—1 



_i ^0 



dtidt 2 <V+(fi)V+(t 2 )> 



x [D^ndti — *a|) — D nm (\t - ti - i 2 |)|l01) 

is the function which controls the Cooperon decay in 
time, i.e. describes electron decoherence for our Id ar- 
ray of quantum dots. The WL correction G\yl m the 
presence of electron-electron interactions is recovered by 
substituting the result ([TOOl into Eq. (1591 . 

Since the behavior of the latter formula was already 
analyzed in details earlier there is no need to repeat this 
analysis here. The dephasing time t v can be extracted 
from the equation J-(t v ) = 1. From Eq. (llOip with a 
good accuracy we obtain 



1 



N ■ 



N-l 
n.m—1 



dr{V+(r)V+(0))D nm (r). (102) 



Combining this formula with Eqs. (1^71) and (|9"9")l , in the 

most interesting limit T — > and for to 3> R(AC + C g ) 
we find 



1 



1 



N-l 



2g TD (N 



— Y In 
1 Z - 



2e J 



9=1 



5 (AC (1 - cos ^)+C g )' 
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which yields 



2gr D 



4tt 



\n(4E c /5) 5\n{4E c /S) 



(103) 



where E c = e 2 /2C g for C g > C and E c = e 2 /4C in the 
opposite case C g <C. 

In order to determine the dephasing length L v = 
\JDt v let us define the diffusion coefficient 

d 2 d 2 gS 



D = 



2tt> 47T 



(104) 



where d = V 1 / 3 is the average dot size. Combining Eqs. 
([TU31 and (ITIHl) . at T = we obtain 



L 



ipO 



= = d^g/\n(4E c /5). 



(105) 



At non-zero T thermal fluctuations provide an addi- 
tional contribution to the dephasing rate 1/t v . Again 
substituting Eqs. (JSTJ) and (jM]) into (1TU2")) . we get 



r v (T) 



1 

T V 



(106) 



where N v = L v /d ~ ^Jt^/td is the number of quantum 
dots within the length L y . We observe that for suffi- 
ciently small N < Ntp (but still N 3> 1) the dephasing 
rate increases linearly both with temperature and with 



the number N. At larger N > \J g / h\[AEc / S\ and/or at 
high enough temperatures N v becomes smaller than N 
and Eq. (|106l) for t v should be resolved self-consistently. 
In this case we obtain 



(3sv / ^D/ 7rT ) 



2/3 



(107) 



thus reproducing the well known result^. Eq. 
(|106p also allows to estimate the temperature T* ~ 
2irg/[r v o min{7V, N v }} at which the crossover to the 
temperature-independent regime (|103[) occurs. We find 



rj-i 



3\n[AE c /S] 
2ttNt d ' 

3m 3/2 [4£ c /<5] 
2nT Dy /g 



N < 



9 



\n[AE c /5] 



N > 



ln[AE c /5 



, (108) 



2. Good metals and granular conductors 

The above analysis and conclusions can be generalized 
further to the case 2d and 3d structures. This gener- 
alization is absolutely straightforward (see, e.g.^) and 
therefore is not elaborated here. At T — » one again 
arrives at the same result for t v0 (|103[) . 

Now we discuss the relation between our present results 
and those derived earlier for weakly disordered metals by 
means of a different approach^. Let us express the dot 



mean level spacing via the average dot size d as <5 = 
l/Nod 3 (where No — mppfa-K 2 is the electron density of 
states at the Fermi level). Then we obtain 



D 



4nN d 



(109) 



Below we consider two different physical limits of (a) 
good metals and (b) strongly disordered (granular) con- 
ductors. For the model (a) we assume that quantum 
dots are in a good contact with each other. In this case g 
scales linearly with the contact area A = jd 2 , where 7 is 
a numerical factor of order (typically smaller than) one 
which particular value depends on geometry. For weakly 
disordered metals most conducting channels in such con- 
tacts can be considered open. Hence, g = p 2 7 A/2n and 



D = jv F d/4, 



(110) 



i.e. D oc d. If most channels are not fully transparent, 
then the factor 7 in (|110j) also accounts for their transmis- 
sions. Comparing Eq. (|110[) with the standard definition 
of D for a bulk diffusive conductor, D = Vi?//3, we im- 
mediately observe that within our model the average dot 
size is comparable to the elastic mean free path, I ~ 7c?, 
as it should be for weakly disordered metals. 
Expressing t v q (|103[) via D, in this limit we get 



•tf>0 



64 m 2 D 3 
^f4ln(D/D c i) 



(111) 



where m is the electron mass and D c \ is constant which 
depends on Ec- Estimating, e.g., Ec ~ e 2 /2d, one ob- 
tains D~i = 4e v / 2lVo/7VF. 

Note that apart from an unimportant numerical pre- 
factor and the logarithm in the denominator of Eq. (|111[) 
the latter result for t v o coincides with that derived for 
a bulk diffusive metal within the framework of a com- 
pletely different approach^. Within that approach lo- 
cal properties of the model remain somewhat ambigu- 
ous and, hence, in the corresponding integrals ini we 
could not avoid using an effective high frequency cutoff 
procedure. This cutoff yields the correct leading depen- 
dence t v q cx D 3 and it only does not allow to recover an 
additional logarithmic dependence on D in (jllip . Our 
present approach is divergence-free and, hence, it does 
not require any cutoffs. 

We can also add that Eq. (|103|) also agrees with our 
earlier results^ derived for quasi-ld and quasi-2d metallic 
conductors. Provided the transversal size a of our array 
is smaller than d one should set A ~ da for 2d and A ~ a 2 
for Id conductors. Then Eq. (|103|) yields t v0 cx D 2 / InZ? 
and t v q oc D/hxD respectively in 2d and Id cases. Up 
to the factor In D these dependencies coincide with ones 
derived previously^. 

Now let us turn to the model (b) of strongly disordered 
and/or granular conductors. In contrast to the situa- 
tion (a), we will assume that the contact between dots 
(grains) is rather poor, and inter-grain electron transport 
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FIG. 11: Ring composed of N quantum dots 



may occur only via limited number of conducting chan- 
nels. In this case the average dimensionless conductance 
g can be approximated by some ,4-independent constant 
g = g c . Substituting g c instead of g into Eq. (I109|) we 
observe that in the case of strongly disordered structures 
one can expect Dal jd. Accordingly, for r^o (I103[) one 
finds 



32n 2 N 2 D 3 ln(D c2 /D) 



(112) 



where D C 2 again depends on Eq- For Eq ~ e 2 /2d we 
have D~2 = 2n\/2N /eg c . Hence, the dependence of t v q 
on D for strongly disordered or granular conductors (|1 12[) 
is it qualitatively different from that for sufficiently clean 
metals ([TTT]) . 

One can also roughly estimate the crossover between 
the regimes (a) and (b) by requiring the values of D = 
^vpd/A (jllOp and D = g c /AirNod to be of the same order. 
This condition yields (ppd) 2 ~ 2irg c /^ 1 and we arrive at 
the estimate for D at the crossover 



D 



o.m 



(113) 



Here we restored the Planck constant h set equal to unity 
elsewhere in our paper. 



3. Ring composed of quantum dots 

Now let us turn to a ring-shaped nanostructure as 
shown in Fig. [TT] For simplicity we will consider the 
case of identical quantum dots (with mean level spacing S 
and dwell time tjj = 2n / {g5)) coupled by junctions with 
conductances gt and the Fano- factor /3 t . Leads are cou- 
pled to the ring at the dots with numbers 1 and L + 1 by 
junctions with conductance g. The interference correc- 
tion to the conductance of n-th junction 5G n was already 
derived in Sec. 2 by means of the non-linear sigma-model 
approach. We obtain 



SG r , 



e 2 g t S 
' 4tt 2 



dt\p t Cn,n+i(t)e 



the aid of the Kirchhoff's law. 
considered here one finds 



For the case Ng <c gt 



5G = 



NL(N - L)g 2 
(2Ng t + L(N - L)g) 



;Sg 



L(N- L)g 2 
*Ng 2 



SG n . 



(115) 

Further procedure is analogous to that implemented 
above for Id arrays. The main difference of the present 
ring-shaped geometry just concerns the form of diffusons 
D mn (t), cooperons Cmn(t) and the fluctuating voltage 
correlators F mn (t) = (V+(t)V+ (0)) v + . We obtain 

N ^ J 2tt -iutd + e[q) 



N 



du) e -»wt+-^ a (m-n) 



c(°) (t) = — V 
mnU N Z^J 27r -iujtd + e(q - 2$/$ ) : 



(117) 



and 



9=1 J 71 



2T uj 2 T 2 D +e 2 (q) ' 



where 



f(g) 



gtrpe 2 e(q) 

n (4Ce(q) + C g ) 2 - 



9tTD& 



e{q) 



IT 



ACe{q) + C g 



(118) 
(119) 

(120) 



cos 2 jp . As above, here C and C g denote 



and e(q) = 1 

respectively the junction and the dot capacitances. 

The above equations are sufficient to evaluate the func- 
tion !F{t) in a general form. Here we are primarily in- 
terested in AB oscillations and, hence, we only need to 
account for the flux-dependent contributions determined 
by the electron trajectories which fully encircle the ring 
at least once. Obviously, one such traverse around the 
ring takes time t > N 2 td- Hence, the behavior of the 
function F(t) only at such time scales needs to be stud- 
ied for our present purposes. In this long time limit F(t) 
is a linear function of time with the corresponding slope 



T'(t > N 2 t d ) 

2e 2 T 2 N-l 



N 



q=l 



did 



f(q)e(q)u coth 



2T 



2tt (w 2 



' D 



i( q ))(uj 2 



2 (q)) 



(121) 



+ (1 - Pt)(C n , n (t) + C n+1 ,n+l{t)) + 



47T7>j- 



(114) 



where C m .„(t) is the Cooperon. The quantum correction 
to conductance of the whole system can be obtained with 



This observation implies that at such time scales electron- 
electron interactions yield exponential decay of the 
Cooperon in time 



C mn (t) » C£l(t)e 



(122) 



where 



T <t> 



= F\t> N 2 t d ) 



(123) 



is the effective dephasing time for our problem. In the 
case C g 3> C and td ~> trc = ^ 7r C g /(e 2 g t ) from Eq. 
(|124p we obtain 



ttNT 
3g t 



T « 1/Ntd, 
T » 1/Nt d , 



(124) 



where Ec = e 2 /(2C g ). These expressions are, of course, 
fully consistent with the results derived above in the case 
of Id chains of quantum dots and weakly disordered dif- 
fusive conductors, cf. also^. 

Let us emphasize again that the above results for T{t) 
apply at sufficiently long times which is appropriate in 
the case of AB conductance oscillations. At the same 
time, other physical quantities, such as, e.g., weak local- 
ization correction to conductance can be determined by 
the function F{t) at shorter time scales. Our general re- 
sults allow to easily recover the corresponding behavior 
as well. For instance, at T td and t <C N 2 tjj we get 



AT 



2tt 
Td 



1/2 



t 3/2 + 



(125) 



in agreement with the results 2 ^. This expression yields 
the well known dependence oc T~ 2 / 3 which - in con- 
trast to Eq. (|124j) - does not depend on N and remains 
applicable in the high temperature limit. 

To proceed further let us integrate the expression for 
the Cooperon over time. We obtain 



C mn {t)dt 



TP \ ^ 



N 



e jv 



'■ (m—n) 



N =i <1 ~ 2$ / $ o) + td/t^ + g/{g t N) 



, (126) 



where the term g/(gtN) in the denominator accounts for 
the effect of external leads and remains applicable as long 
as Ng < g t . Combining Eqs. (TiTI)) , ([TT5j) and (Tl2"6")) after 
summation over q we arrive at the final result 



§g ab = e 2 L(N-L)g 2 



2nNg 2 



-N 



cos(47T$/$ )) 



Va 2 - l(z N + z- N - 2 cos(4tt$/$o)) ' 



(127) 



where a = 1 - 



t> + gTJv anc ^ z = a + V« 2 — I- This equa- 
tion with Eq. (1124)) fully determines AB oscillations of 
conductance in nanorings composed of metallic quantum 
dots in the presence of electron-electron interactions. 
Expanding Eq. (|127l) in Fourier series we obtain 



5G AB = SG(k) cos (47r/fc$/$ ) 



(128) 
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FIG. 12: Temperature dependence of the first three harmonics 
of AB conductance for gt = 500, g = 30, N = 10, /3t = 1 and 
td/trc = 120. 



where 



SG (k) = e 2 L(N-L)g 2 (p t a + l-!3 t ) ^ Nlkl 
2irNg 2 v / a 2 - 1 



(129) 



In the limit ^> tjj we have z « 1 + \l 2td / t$ 
hence SG^ behaves as 



6G (k) 



-N\k\. 



oc e 



(130) 



i.e. at hight temperatures log \SG\ scales with N as N 3 / 2 
while at low temperatures it scales as N. The tempera- 
ture dependence of the first three harmonics of AB con- 
ductance in the presence of electron-electron interactions 
is depicted in Fig. 12. 

The results obtained here allow to formulate quanti- 
tative predictions regading the effect of electron-electron 
interactions on Aharonov-Bohm oscillations of conduc- 
tance for a wide class of disordered nanorings embraced 
by our model. Of particular interest is the situation of 
large number of dots N ^> 1 which essentially mimics the 
behavior of diffusive nanostructures. In order to establish 
a direct relation to this important case it is instructive 
to introduce the diffusion coefficient D = d 2 /(2Tu) and 
define the electron density of states v = l/(d 3 S), where 
d is a linear dot size. Then we obtain with exponential 
accuracy: 



5G^ 



e -\k\(C/C^) 
e -|fe|(£/£^) 3 / 2 



T< D/(Cd), 
T^>D/(Cd). 



Here we introduced the ring perimeter C = Nd and the 
effective decoherence length 



— ro^ 



In 



1/2 



1/3 



T < D/(£d), 
T > D/{Cd). 



k=l 



Note in the high temperature limit T D / (Cd) the 
above results match with those derived earlier for metal- 
lic nanorings with the aid of different approache o 39 ' 40 . 
On the other hand, at lower T our results are different. 
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This difference is due to low temperature saturation of 
which was not accounted for i n 39 ' . A non-trivial feature 
predicted here is that - in contrast to weak localization 4 
- the crossover from thermal to quantum dephasing is 
controlled by the ring perimeter C. This is because only 
sufficiently long electron paths fully encircling the ring 
are sensitive to the magnetic flux and may contribute to 
AB oscillations of conductance. 

We believe that the quantum dot rings considered here 
can be directly used for further experimental investiga- 
tions of quantum coherence of interacting electrons in 
nanoscale conductors at low temperatures. 

IV. COMPARISON WITH EXPERIMENTS AND 
CONCLUDING REMARKS 

Turning to the experimental situation in the field, it is 
important to emphasize again that low temperature sat- 
uration of the electron decoherence time has been repeat- 
edly observed in numerous experiments and is presently 
considered as firmly established and indisputably existing 
phenomenon. Although in some cases this phenomenon 
can be attributed to various extrinsic mechanisms, like 
magnetic impurities, overheating etc., in the vast major- 
ity of cases none of such extrinsic mechanisms can rea- 
sonably account for experimental observations. On the 
other hand, it was demonstrated above that electron- 
electron interactions universally provide non-vanishing 
electron dephasing down to T = in all types of dis- 
ordered conductors. Therefore, it would be interesting 
to perform quantitative comparison between our univer- 
sal formula for t v q, Eq. (|103[) . and experimental values 
of the electron decoherence time measured in different 
structures. 

Note that in some of our earlier publications^ ' 41 ' 42 we 
have already demonstrated a good quantitative agree- 
ment between our theoretical predictions^ and experi- 
mental data for t v q obtained for numerous metallic wires 
and quasi-ld semiconductors. Here we address the ex- 
periments on quantum dot structures as well as on both 
weakly and highly disordered metals. 

First turning to quantum dots, we recall that in all 
14 samples reported in experiments with open quantum 
dots performed by different groups^ 3 - - — the values t v0 
were found to rather closely follow a simple dependence^, 

7^0 «td. (131) 

This approximate scaling was observed within the inter- 
val of dwell times td of about 3 decades, see Fig. 5 
in 4 ^. Our Eq. ()103[) essentially reproduces this scaling, 
especially having in mind that the dimensionless conduc- 
tance g was of order one (or slighlty larger) in almost all 
samples^^. To the best of our knowledge no alter- 
native explanation for the scaling (|131l) has been offered 
until now. Thus, we conclude that our theory is clearly 
consistent with the available experimental data on zero 
temperature electron dephasing in open quantum dots. 
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FIG. 13: The low temperature dephasing times observed in 
various experiments for the following samples: Au-1 to Au- 
6 6 , Au-7 50 , Au-8 and Au-lC 51 (■); 44 samples (AuPd and 
AgPd) 47 (o); 18 samples 48 : Au 2 Al (<*), Sb (<), Sc 85 Ag 15 (►), 
V 3 A1 (B); 9 samples (CuGeAu)^ (ffl) ; 15 samples (Au, Ag 
and Cu)S an d AgMI6N 53 (0); CF-1 and CF-2 54 (A); A, B 
(Aup-, Aul- 7 -, Agl^ and Ag3^ (□); S, M and L (Pt)& (.); 
D and F 59 (>); Ag, AgFel and AgFe3^ (V); 10 sampled 
within the box (E2); 2 (Au) e2 and 1 (Au)& (a); al-1^ (♦). 
Our Eq. (Illip for 7 = 0.2 and 1 is indicated respectively 
by solid and dashed lines, while Eq. (|112|) for g c = 150 is 
depicted by dashed-dotted line. 



Let us now consider spatially extended disordered con- 
ductors. As our theory of dephasing by electron-electron 
interactions predicts a rather steep increase of r v o with 
the system diffusion coefficient D, for most weakly dis- 
ordered metals as t v q oc D 3 , we can conclude that for a 
large number of disordered conductors t v q strongly in- 
creases with increasing D. This trend is indeed quite 
obvious for relatively weakly disordered conductors. On 
the other hand, Lin and coworker o 9 ' 47 " — analyzed numer- 
ous experimental data for t v q obtained by various groups 
in rather strongly disordered conductors with D < 10 
cm 2 /s and observed systematic decrease of t v o with in- 
creasing D. The data could be fitted by the dependence 
t v q D~ a with the power a > 1. This trend is clearly 
just the opposite to one observed in less disordered con- 
ductors with D > 10 cm 2 /s. 

In Fig. 13 we collected experimental data for t v q 
obtained in about 130 metallic samples with similar 
Fermi velocities and diffusion coefficients varying by ~ 3 
decades, from D 0.3 cm 2 /s to D ss 350 cm 2 /s. The 
data were taken from about 30 different publications 
listed in the figure caption. We see that the the measured 
values of t v q strongly depend on D. Furthermore, this 
dependence turns out to be non-monotonous: For rela- 
tively weakly disordered structures with D > 10 cm 2 /s 
t v o increases with increasing D, while for strongly disor- 
dered conductors with D < 10 cm 2 /s the opposite trend 
takes place. In addition to the data points in Fig. 13 
we indicate the dependencies t v q(D) (|llip and (I112[) for 
two models (a) and (b) discussed above. 
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We observe that for D > 10 cm 2 /s the data points 
clearly follow the scaling (|lllj) . Practically all data 
points remain within the strip between the two lines cor- 
responding to Eq. (jlllj) with 7 = 1 (dashed line) and 
7 = 0.2 (solid line). On the other hand, for more disor- 
dered conductors with D < 10 cm 2 /s the data are consis- 
tent with the scaling (|112[) obtained within the model (b) . 
We would like to emphasize that theoretical curves (|111[) 
and (|112|) are presented in Fig. 13 with no fit parameters 
except for a geometry factor 7 for the first dependence 
and the value g c ps 150 for the second one. This value of 
g c was estimated from the crossover condition (j!13p with 
D ~ 10 cm 2 /s and 7 ~ 1. 

Now let us consider the data for strongly disordered 
conductors with D < 10 cm 2 /s. As we already pointed 
out, the agreement between the data and the dependence 
(|112p predicted within our simple model (b) is reasonable, 
in particular for samples with D < 3 cm 2 /s. At higher 
diffusion coefficients most of the data points indicate a 
weaker dependence of r^o on D which appears natural 
in the vicinity of the crossover to the dependence (|llll) . 
The best fit for the whole range 0.3 cm 2 /s < D < 10 
cm 2 /s is achieved with the function t v q oc D~ a with the 
power a ~ 1.5 -7- 2. 

Thus, we conclude that our theory allows to qualita- 
tively understand and explain seemingly contradicting 
dependencies of t v o on D observed in weakly and strongly 
disordered conductors. While the trend "less disorder - 
less decoherence" (|111[) for sufficiently clean conductors 
is quite obvious, the opposite trend "more disorder - less 
decoherence" in strongly disordered structures requires a 
comment. The latter dependence may indicate that with 
increasing disorder electrons spend more time in the ar- 
eas with fluctuating in time but spatially uniform poten- 
tials. As we already discussed in the beginning of Sec. 3, 
such fluctuating potentials do not dephase and thus t v0 
gets effectively increased. In other words, in this case 
the corresponding dwell time td m Eq. (|103|) becomes 
longer with increasing disorder and, hence, the electron 
decoherence time t v q does so too. 

Note that since local conductance fluctuations increase 



with increasing disorder, several grains can form a cluster 
with internal inter-grain conductances strongly exceeding 
those at its edges. In this case fluctuating potentials re- 
main almost uniform inside the whole cluster which will 
then play a role of an effective (bigger) grain/dot. Ac- 
cordingly, the average volume of such "composite dots" 
V oc 1/5 may grow with increasing disorder, electrons 
will spend more time in these bigger dots and, hence, the 
electron decoherence time (|103[) will increase. 

The above comparison with experiments confirms that 
our previous quasiclassical results 4 for t v o are applicable 
to relatively weakly disordered structures with D > 10 
cm 2 /s, while for conductors with stronger disorder differ- 
ent expressions for t v0 (e.g., Eq. (|112[) ) should be used. 
Our analysis also allows to rule out scattering on mag- 
netic impurities as a cause of low temperature saturation 
of t v . The latter mechanism can explain neither strong 
and non-trivial dependence of the electron decoherence 
time on D nor even the level of dephasing observed in 
numerous experiments. E.g., in order to be able to at- 
tribute dephasing times as short as t v q < 10~ 12 s to 
magnetic impurities one needs to assume huge concen- 
tration of such impurities ranging from few hundreds to 
few thousands ppm which appears highly unrealistic, in 
particular for systems like carbon nanotubes, 2DEGs or 
quantum dots. Similar arguments were independently 
put forward by Lin and coworker a 47 ' 49 . 

Thus, although electron dephasing due to scattering 
on magnetic impurities is by itself an interesting issue, 
its role in low temperature saturation of t v in disordered 
conductors is sometimes strongly overemphasized. Since 
the latter phenomenon has been repeatedly observed in 
all types of disordered conductors, the physics behind it 
should most likely be universal and fundamental. We 
believe - and have demonstrated here - that it is indeed 
the case: Zero temperature electron decoherence in all 
types of conductors discussed above is caused by electron- 
electron interactions. 

This work was supported in part by RFBR grant 09- 
02-00886. 
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